HAL 

archives-ouvertes r 


Probabilistic initial value problem for cellular 

automaton rule 172 

Henryk Fuks 


► To cite this version: 

Henryk Fuks. Probabilistic initial value problem for cellular automaton rule 172. Automata 2010 - 
16th Inti. Workshop on CA and DCS, 2010, Nancy, France, pp.29-40. hal-01185497 


HAL Id: hal-01185497 
https://hal.inria.fr/hal-01185497 

Submitted on 20 Aug 2015 


HAL is a multi-disciplinary open access 
archive for the deposit and dissemination of sci¬ 
entific research documents, whether they are pub¬ 
lished or not. The documents may come from 
teaching and research institutions in France or 
abroad, or from public or private research centers. 


L ’archive ouverte pluridisciplinaire HAL, est 
destinee au depot et a la diffusion de documents 
scientifiques de niveau recherche, publies ou non, 
emanant des etablissements d’enseignement et de 
recherche frangais ou etrangers, des laboratoires 
publics ou prives. 



Automata 2010 — 16'^ Inti. Workshop on CA and DCS 


DMTCS proc. AL, 2010, 29^0 


Probabilistic initial value probiem for ceiluiar 
automaton rule 172 


Henryk Fuks ^ 

Department of Mathematics, Brock University, St. Catharines, ON L2S 3A1, Canada 


We present a method of solving of the probabilistic initial value problem for cellular automata (CA) using CA rule 
172 as an example. For a disordered initial condition on an infinite lattice, we derive exact expressions for the 
density of ones at arbitrary time step. In order to do this, we analyze topological structure of preimage trees of 
finite strings of length 3. Level sets of these trees can be enumerated directly using classical combinatorial methods, 
yielding expressions for the number of n-step preimages of all strings of length 3, and, subsequently, probabilities of 
occurrence of these strings in a configuration obtained from the initial one after n iterations of rule 172. The density 
of ones can be expressed in terms of Fibonacci numbers, while expressions for probabilities of other strings involve 
Lucas numbers. Applicability of this method to other CA rules is briefly discussed. 

Keywords: cellular automata, initial value problem, preimage trees 


1 Introduction 

While working on a certain problem in complexity engineering, that is, trying to construct a cellular au¬ 
tomaton rule performing some useful computational task, the author encountered the following question. 
Let / : {0,1}^ —1- {0,1} be defined as 



( 1 ) 


This function may be called selective copier, since it returns (copies) one of its inputs X 2 or X 3 depending 
on the state of the first input variable xi. Suppose now that s be a bi-infinite sequence of binary symbols, 
i.e., s = ... S- 2 S- 1 S 0 S 1 S 2 ..., i C K. We will transform this string using the selective copier, that is, 
for each i, we keep Si if it is preceded by 0, or replace it by otherwise, so that each Si is simulta¬ 
neously replaced by f{si-i, Si, Si+i). Consider now the question: Assuming that the initial sequence is 
randomly generated, what is the proportion of 1 ’s in the sequence after n iterations of the aforementioned 
procedure? 

iThe author acknowledges partial financial support from Natural Sciences and Engineering Research Counclil of Canada, in the 
form of a Discovery Grant. 
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Function defined by eq. Q is a local function of cellular automaton rule 172, using Wolfram number¬ 
ing, and the aforementioned question is an example of a broader class of problems, which could be called 
probabilistic initial value problems for cellular automata: given initial distribution of infinite configura¬ 
tions, what is the probability of occurrence of a given finite string in a conhguration obtained from the 
initial one by n iterations of the cellular automaton rule? In what follows, we will demonstrate how one 
can approach probabilistic initial value problem using cellular automaton rule 172 as an example. 

2 Basic definitions 

Let Q — {0,1, ...TV— 1} be called a symbol set, and let S{Q)he the set of all bisequences over Q, where by 
a bisequence we mean a function on K to Q. Set iS(f/) will be called the configuration space. Throughout 
the remainder of this text we shall assume that Q = {0,1}, and the configuration space S{Q) = {0,1}^ 
will be simply denoted by S. 

A block of length n is an ordered set bobi... b„-i, where n G IM, bi G G- Let n S N and let Bn denote 
the set of all blocks of length n over G and B be the set of all finite blocks over G- 

For r G N, a mapping / : {0, i—>. {Q, 1} will be called a cellular automaton rule of radius r. 

Alternatively, the function / can be considered as a mapping of B 2 r+i into Bq = G = {0,^}- 

Corresponding to / (also called a local mapping) we dehne a global mapping F : S —t S such that 
(F(s))j = /(si-r,..., Si,..., Si+,.) for any s gS. 

A block evolution operator corresponding to / is a mapping i ■. B B defined as follows. Let r G N 
be the radius of /, and let a = oqUi ... On-i G Bn where n > 2r -|- 1 > 0. Then 

f(n) = Oi+l: ■ ■ ■ I (2) 

Note that if 6 G i? 2 r+i then f{b) = f(6). 

We will consider the case of G = {0,1} and r = 1 rules, i.e., elementary cellular automata. In this 
case, when b G B 3 , then f{b) = f (5). The set B 3 = {000,001,010, 011,100,101,101,110, 111} will be 
called the set of basic blocks. 

The number of n-step preimages of the block b under the rule / is dehned as the number of elements 
of the set f“"(6). Given an elementary rule /, we will be especially interested in the number of n-step 
preimages of basic blocks under the rule /. 

3 Probabilistic initial value problem 

The appropriate mathematical description of an initial distribution of conhgurations is a probability mea¬ 
sure p, on S. Such a measure can be formally constructed as follows. If 6 is a block of length k, i.e., 
b = bgbi ... bk-i, then for i gTLsno dehne a cylinder set. The cylinder set is a set of all possible conhg¬ 
urations with hxed values at a hnite number of sites. Intuitively, measure of the cylinder set given by the 
block b = bo ... bk-i, denoted by p[Ci{b)], is simply a probability of occurrence of the block h in a place 
starting at i. If the measure p is shift-invariant, than p(Ci{b)) is independent of i, and we will therefore 
drop the index i and simply write p{C{b)). 

The Kolmogorov consistency theorem states that every probability measure p satisfying the consistency 
condition 


p[Cfibi... bk)] = p[Cfibi ...bk, 0)] -b p[Ci{bi ...bk, 1)] 


(3) 
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extends to a shift invariant measure on S (Dynkin 1969|l .For p G [0,1], the Bernoulli measure dehned 


as p,p[C{b)] = p^{l — p)^~f where j is a number of ones in b and k — j is a number of zeros in b, is 
an example of such a shift-invariant (or spatially homogeneous) measure. It describes a set of random 
conhgurations with the probability p that a given site is in state 1. 

Since a cellular automaton rule with global function F maps a configuration in S to another configura¬ 
tion in S, we can dehne the action of F on measures on S. For all measurable subsets of 5 we dehne 
{Fp){E) = p.{F~^{E)), where F~^{E) is an inverse image of E under E. 

If the initial configuration was specihed by p,p, what can be said about F'^p.p (i.e., what is the probability 
measure after n iterations of F)1 In particular, given a block b, what is the probability of the occurrence 
of this block in a conhguration obtained from a random conhguration after n iterations of a given rule? 

The general question of hnding the iterrates of the Bernoulli measure under a given CA has been 


extensively studied in recent years by many authors, including, among others. 

Lind 

00 

Ferrari et al. 

(2000|l;|Maass and Martinez ( 

2003)1; |Host et al.|(|20031l;|Pivato and Yassawi|( 

2002 

20041 

Maass et al. 


and 


Maass et al. (2006|l. In this paper, we will approach the problem from somewhat different 


angle, using very elementary methods and without resorting to advanced apparatus of ergodic theory and 
symbolic dynamics. We will consider iterates of the Bernoulli measure by analyzing patterns in preimage 
sets. 

For a given block b, the set of n-step preimages is f“"(6). Then, by the definition of the action of F on 
the initial measure, we have 

{F-pp){C{b))=pp{F--{Cib))), (4) 

and consequently 

(F>p)(C(&))= ^ pp{a). (5) 

Let us dehne the probability of occurrence of block 5 in a conhguration obtained from the initial one by n 
iterations of the CA rule as 

Pn{b) = (F"pp)(C(6)). (6) 

Using this notation, eq. 0 becomes 


Pn{b) = 




Po{a). 


(7) 


If the initial measure is pi/ 2 , then all blocks of a given length are equally probable, and Po{a) = 
where |a| is the length of the block a. For elementary CA rule, the length of n-step preimage of b is 
2 n + |6|, therefore 

P„(6) = 2-l^l-2"cardf-”(6). (8) 

This equation tells us that if the initial measure is symmetric (p. 1 / 2 ), then all we need to know in order to 
compute Pn{b) is the cardinality of f“”(6). One way to think about this is to draw a preimage tree for b. 
We start form & as a root of the tree, and determine all its preimages. Then each of these preimages is 
connected with b by an edge. They constitute level 1 of the preimage tree. Then, for each block of level 1, 
we again compute its preimages and we link them with that block, thus obtaining level 2. Repeating this 
operation ad infinitum, we obtain a tree such as the one shown in Figure [T] In that hgure, hve levels of the 
preimage tree for rule 172 rooted at 101 are shown, with only first level labelled. 
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Fig. 1: Preimage tree for rule 172 rooted at 101. 

Note that cardf“"(6) corresponds to the number of vertices in the n-th level of the preimage tree. 
One thus only needs to know cardinalities of level sets in order to use eq. while the exact topology 
of connections between vertices of the preimage tree is unimportant. The key problem, therefore, is to 
enumerate level sets. In order to answer the question posed in the introduction, we need to compute Pn (1) 
for rule 172, which, in turn, requires that we enumerate level sets of a preimage tree rooted at 1. It turns 
out that for rule 172 the preimage tree rooted at 1 is rather complicated, and that it is more convenient to 
consider preimage trees rooted at other blocks. In the next section, we will show how to express Pn(l) by 
some other block probabilities. From now on, f will exclusively denote the block evolution operator for 
rule 172. 

4 Block probabilities 

Since f-i(l) = {010,011,101,111}, we have P„+i(l) = P„(010) + P„(011) + P„(101) + P„(lll). 
Due to consistency conditions (eq. [^, P„(010) + P„(011) = P„(01), and we obtain 


Pn+l(l) = Pn(Ol) + P„(101) + P„(lll). 


(9) 
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P„(l) = P„_i(001) + 2P„_i(101) + P„_i(lll). (10) 

By using eq. (j^ and defining c„ = T’n(l) we obtain 


card f-”+i (001) + 2cardf-”+i(101) + cardf-”+i(lll) 

22n+l 


( 11 ) 


This means that in order to compute c„, we need to know cardinalities of n-step preimages of 001, 101, 
and 111. 


5 Structure of preimage sets 

The structure of level sets of preimage trees rooted at 001, 101, and 111 will be described in the following 
three propositions. 

Proposition 5.1 Block b belongs to f“”(001) if and only if it has the structure 

6 = ★ ★ . 001 ^ 'k . . , ( 12 ) 


where * represents arbitrary symbol from the ief {0,1}. 

Let us first observe that f“^(001) = {00010,00011,10010,10011}, which means that f~^(001) can be 
represented as **001 **. Similarly, therefore, f~^(001) has the structure ***001***, and by induction, 
for any n, the structure of f“”(001) must be * * ^.. * 001 ** ■ ■ -f - LI 

71 71 

Proposition 5.2 Block b belongs to f“”(101) if and only if it has the structure 

b = * * ^.. * ■ ■ ■ ay^llOl, (13) 

n—1 

where ai £ {0, l}/or i = 1 ,... ,n and the string aia 2 ■ ■ ■ an does not contain any pair of adjacent zeros, 
that is. OiOi+i f 00 for all i = 1 ,... ,n — 1. 

Two observations will be crucial for the proof First of all, f~^(101) = (01101,11101}, thus f~^(101) 
has the structure *1101. Furthermore, we have f~^(1101) = (011101,101101,111101}, meaning that if 
1101 appears in a configuration, and is not preceded by 00, then after application of the rule 172, 1101 
will still appear, but shifted one position to the left. All this means that if b is to be an n-step preimage of 
101, it must end with 1101. After each application of rule 172 to b, the block 1101 will remain at the end 
as long as it is not preceded by two zeros. 

Now, let us note that f“^(00) = (0000,0001,1000,1001,1100}, which means that preimage of 00 is 
either 1100 or *00*. Therefore, we can say that if 00 is not present in the string aia 2 ... a„, it will not 
appear in its consecutive images under f. Thus, block 1101 will, after each iteration of /, remain at the 
end, and will never be preceded by two zeros. Eventually, after n iterations, it will produce 101, as shown 
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in the example below. 


10 10 0 1 
110 0 1 
0 0 0 1 
0 0 1 

0 1 

1 


110 110 1 
10 110 1 
0 110 1 
110 1 
1 0 1 
0 1 


What is left to show is that not having 00 in ai 02 ... is necessary. This is a consequence of the fact 
that f(*00*) = 00, which means that if 00 appears in a string, then it stays in the same position after 
the rule 172 is applied. Indeed, if we had a pair of adjacent zeros in 0102 ... a„, it would stay in the 
same position when f is applied, and sooner or later block 1101, which is moving to the left, would come 
to the position immediately following this pair, and would be destroyed in the next iteration, thus never 
producing 101. Such a process is illustrated below, where after three iterations the block 1101 is destroyed 
due to “collision” with 00. □ 


10 10 0 
110 0 
0 0 0 
0 0 
0 


0 110 110 1 
0 10 110 1 
0 1110 1 
0 110 1 

0 10 1 
0 1 1 


Proposition 5.3 Block b belongs to f ”(111) if and only if it has the structure 


b — ^ -k . * 0102 . . . On+S, 
n—2 


(14) 


where at G {0, l}/or i = 1 ,... ,n and the string 0102 ... a„ satisfies the following three conditions: 

(i) aiOi^i f 00 for all i = 2 .. .n + 4; 

(ii) f 110 U^_|_2G.n+3ttn+4 7^ 110, 

(Hi) if aia 2 f 00, then an+ia„+ 2 a „+3 f 110. 

We will present only the main idea of the proof here, omitting some tedious details. It will be helpful to 
inspect spatiotemporal pattern generated by rule 172 first, as shown in Figure]^ Careful inspection of this 
pattern reveals three facts, each of them easily provable in a rigorous way: 

(FI) A cluster of two or more zeros keeps its right boundary in the same place for ever. 

(F2) A cluster of two or more zeros extends its left boundary to the left one unit per time step as long as 
the left boundary is preceded by two or more ones. If the left boundary it is preceded by 01, it stays 
in the same place. 


(F3) Isolated zero moves to the left one step at a time as long as it has at least two ones on the left. If an 
isolated zero is preceded by 10, it disappears in the next time step. 
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i 



Fig. 2: Example of a spatiotemporal pattern produced by rule 172. 

Let us first prove that (i)-(iii) are necessary. Condition (i) is needed because if we had 00 in the string 
02 ... a„+ 5 , its left boundary would grow to the left and after n iterations it would reach sites in which 
we expect to hnd the resulting string 111. 

Moreover, string 0102 ... a „+5 cannot have Oil at the end position, one site before the end, or two sites 
before the end. If it had, 0 preceded by two 1 ’s would move to the left and, after n iterations, it would 
reach sites where we want to find 111. The only exception to this is the case when aoOi = 00. In this 
case, even if Oil is in the second position from the end, it will disappear in step n — 1. This demonstrates 
that (ii) and (iii) are necessary. 

In order to prove sufficiency of (i)-(iii), let us suppose that the string b satisfies all these conditions 
yet f"(&) f 111. This would imply that at least one of the symbols of f"(6) is equal to zero. However, 
according to what we stated in F1-F3, zero can appear in a later conhguration only as a result of growth 
of an initial cluster of two of more zeros, or by moving to the left if it is preceded by two ones. This, 
however, is impossible due to conditions (i)-(iii). □. 

6 Enumeration of preimage strings 

Once we know the structure of preimage sets, we can enumerate them. For this, the following lemma will 
be useful. 

Lemma 6.1 The number of binary strings 0102 ... a„ such that 00 does not appear as two consecutive 
terms OiOi+i is equal to Fn+ 2 > where F„ is the n-th Fibonacci number. 

This result will be derived using classical transfer-matrix method. Let g{n) be the number of binary 
strings 0102 ... a„ such that 00 does not appear as two consecutive terms OiOi+i. We can think of such 
string as a walk of length n on a graph with vertices ui = 0 and U 2 = 1 which has adjacency matrix A 
given by All = 0, A 12 = ^21 = A 22 = 1. One can prove that the generating function for g, 

00 

G(A) = ^p(n+l)A", 

n=0 


(15) 
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can be expressed by G'(A) = Gii(A) + Gi 2 (A) + G 2 i(A) + G 22 (A), where 

(—det(J — XA : j, i) 




det(/ — XA) 


(16) 


and where {M : j, i) denotes the matrix obtained by removing the j — th row and i — th column of M. 
Proof of this statement can be found, for example, in |Stanle^ ( |1986| l. Applying this to the problem at hand 
we obtain 

G(A)= 


-1 + A + A2- 

By decomposing the above generating function into simple fractions we get 

igVO 9 , 9 inVO 


G(A) = 


10 


A + V' 


x + i-^p ’ 


where ijj = i + ^v/bis the golden ratio. Now, by using the fact that 

n+1 


and by using a similar expression for 


^ = (zi 

A + V” 

, we obtain 


\n 




(17) 


(18) 


(19) 


G(A) = ^ K+sA” 

n=0 


where Fn is the n-th Fibonacci number, F„ = 


(1-V')” 

TH 


Proposition 6.1 The cardinalities of preimage sets o/OOl, 100, 101 and 111 are given by 

cardf-”(001) = 4", 

cardf-’"(101) = 2 ”- 1 f„+ 2 , 

cardf-"(lll) = 2 "F„+ 3 . 


( 20 ) 


. This implies that g{n) = 7^+2 ■ ^ 


( 21 ) 

( 22 ) 

(23) 


Proof of the first of these formulae is a straightforward consequence of Proposition o We have 2n 
arbitrary binary symbols in the string b, thus the number of such strings must be 2^" = 4". 


The second formula can be immediately obtained using Lemma 6.1 and Proposition 5.1 Since the 


first n — 1 symbols of f "(101) are arbitrary, and the remaining symbols form a sequence of n symbols 
without 00, we obtain 

cardf-"(101) = 2"-ii^„+2. (24) 

In order to prove the third formula, we will use Proposition |5.3| We need to compute the number 
of binary strings 0102 ... a „+5 satisfying conditions (i)-(iii) of Proposition |5. 3 [ Le us first introduce a 
symbol Q!ia 2 • • ■ ctk to denote the string of length k in which no pair 00 appears. Then we define; 


A is the set of all strings having the form q;iQ !2 ... an+ 5 , 
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• Ai is the set of all strings having the form aia 2 ■ ■ ■ a„+ 2 ll 0 , 

• A 2 is the set of all strings having the form aia 2 ■ ■ ■ an+iHOl, 

• A^ is the set of all strings having the form aia 2 ■ ■ ■ a„11010, 

• A 4 is the set of all strings having the form aia 2 ■ ■ ■ a„11011, 

• B is the set of all strings having the form 001aia2 ... an+ 2 , 

• Bi is the set of all strings having the form 001 aia 2 ... a„_illO, 

• i ?2 is the set of all strings having the form 001 q;iq ;2 ... a„_ 2 ll 01 . 


The set fl of binary strings 0402 ■ ■ • 0^+5 satisfying conditions (i)-(iii) of Proposition 5.3 can be now 
written as 

n = A\{AiU A2U A3U A4)U B\{BiU B2). (25) 

Since Ai ... A 4 are mutually disjoint, and Bi and B 2 are disjoint too, the number elements in the set 12 is 


card ft = card A — card Ai — card A 2 — card A 3 — card A 4 

+ card B — card Bi — card B 2 , 


(26) 


which, using Lemma | 6 T[ yields 

card n = Fn+7 — (T)i+4 + Fn+3 + Fn 4-2 + T’n+ 2 ) + Fn+4 ~ (^n+l + Fn). (27) 

Using basic properies of Fibonacci numbers, the above simplihes to card 12 = AFn+ 3 . Now, since in the 
Proposition |5 .3 1 the string oi ... a „+5 is preceded by n — 2 arbitrary symbols, we obtain 


cardf-"(lll) = . 4 F „+3 = 2 ”F„+ 3 , 


(28) 


what was to be shown. 

7 Density of ones 


Using results of the previous section, eq. (Ill can now be rewritten as 

4”-i + 2"-iU„+i + 2"-1F„+2 


Cn — 


which simplifies to 

or, more explicitly, to 


Cn — 


22n+l 

1 Fn+3 


Cn — 


8 2^+2’ 

1 ^ (1 + 75) "+^ - (1 - 75)"+^ 

22n+5v/5 


(29) 


(30) 


(31) 


Wolfram 


(1994 1 . We can 


Obviously, lim„_>oo Cn = g, in agreement with the numerical value reported in 
see that c„ converges toward Coo exponentially fast, with some damped oscillations superimposed over 
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0.835 

- 1 - 1 -^^ -1-1-1-^-1 


numerical 


theoretical 

0.83 


^ 0.825 

- 

o 


^ 0.82 

- 



0.815 

- 

C 


0.81 

. • . • ■ 

• 

0.805 

- 

0.8 

- * -1- ^^ -1-1-1- ^ - 


1 23456789 10 

n 


Fig. 3 : Plot of the ratio as a function of time step n. Numerical results were obtained by iterating rule 

C-n Coo 

172 on a a configuration of length 10® with periodic boundary conditions. 


the exponential decay. This is illustrated in Figure where, in order to emphasize the aforementioned 
oscillations, instead of c„ we plotted the ratio 




Cn+1 Coo 
Cn Coo 


(32) 


as a function of n. One can show that converges to the half of ratio divina (golden ratio), -0/2 « 
0.809016 ..., as illustrated in Figure]^ We can see from this figure that the convergence is very fast and 
that the agreement between numerical simulations and the theoretical formula is nearly perfect. 


8 Further results 


Results obtained in the previous two sections suffice to compute block probabilities for all blocks of 
length up to 3. Proposition |6. 1 [ together with eq. ([^ yields formulas for P„(001), P„(101), and P„(lll). 
Consistency conditions give P„(01) = Pn(OOl) + Pn(lOl). Furthermore Pji(lO) = P„(01) due to the 
fact that P„(10) + P„(00) = P„(01) + P„(00) = P„(0). Applying consistency conditions again we have 
Pn{l) = Pn(lO) +P„(11), hence P„(ll) = P„(l) - P„(10), and, similarly, P„(00) = P„(0) - P„(10). 
This gives us probabilities of all blocks of length 2. Probabilities of blocks of length 3 can be obtained in 
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P„(000) = P„(00)-P„(100), 

P„(110) = P„(11)-P„(111), 

P„(011) = P„(11)-P„(111), 

P„(010) = P„(01) - P„(011). 

The only missing probability, P„(100) is the same as P„(001) because Pn(lOO) + P„(000) = P„(001) + 
P„(000) = P„(00). The following formulas summarize these results. 

P„(000) = 5/8 - 2-’^-2p„+3 - 2-"-V„+2, 

P„(001) = 1/8, 

P„(010) = l/8-2-’^-3p„+i, 

P„(011) = 2-’^-4p„+2, 

P„(100) = 1/8, 

P„(101) = 2-"-4p„+2, 

P„(110) = 2-"-4p„+2, 

P„(lll) = 2-’^-3p„+3, 

where P„ = 2P„+i — P„ is the n-th Lucas number. We can also rewrite these formulas in terms of 
cardinalities of preimage sets using eq. (j^, as stated below. 

Theorem 8.1 Let f be the block evolution operator for CA rule 172. Then for any positive integer n we 
have 

card f“" (000) = 

card f“" (001) = 

card f“" (010) = 

cardf“"(011) = 

card f“"( 100) = 

cardf“"(101) = 

cardf“"(110) = 

cardf“"(lll) = 

where P„ is the n-th Fibonacci number, Fn 
number, L„ = r/" + (1 — r/)”. 

9 Concluding remarks 

The method for computing block probabilities in cellular automata described in this paper is certainly 
not applicable to arbitrary CA rule. It will work only if the structure of level sets of preimage trees is 


5.4"_2"+ip„+3-2"-ip„+2, 


4" - 2"P 


n+l5 


nn—l T 

4", 

2"-ip, 


n+25 
'n+2, 

2"Pn+3, 


nn-1 T 
Z Ijri 


- (1 - r/>)” , , , r- 

- P-, "0 = 4 + 5 V 5, and P„ is the n-th Lucas 

V5 
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sufficiently regular so that the level sets can be enumerated by some known combinatorial technique. 
Altough “chaotic” rules like rule 18, or complex rules such as rule 110 certainly do not belong to this 
category, in surprisingly many cases significant regularities can be detected in preimage trees. Usually, this 
applies to “simple” rules, those which in Wolfram classification belong to class I, class II, and sometimes 
class III. Rule 172 reported here is one of the most interesting among such rules, primarily because the 
density of ones does not converge exponentially to some fixed value as in many other cases, but exhibits 
subtle damped oscillations on top of the exponential decay. Furthermore, the appearance of Fibonacci and 
Lucas numbers in formulas for block probabilities is rather surprising. 

One should add at this point that the convergence toward the steady state can be slower than exponential 


even in fairly “simple” cellular automata. Using similar method as in this paper, it has been found in Fuks 


and Haroutunian (2009 |l that in rule 14 the density of ones converges toward its limit value approximately 


as a power law. The exact formula for the density of ones in rule 14 involves Catalan numbers, and the 
structure of level sets is quite different than the one reported here. Rule 142 exhibits somewhat similar 
behavior too, as reported in Fuks ( 2006| l. 

As a final remark, let us add that the results presented here assume initial measure ^. 1 / 2 - This can 
be generalized to arbitrary pp. In order to do this, one needs, instead of straightforward counting of 
preimages, to perform direct computation of their probabilities using methods based on Markov chain 
theory. Work on this problem is ongoing and will be reported elsewhere. 
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